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ABSTRACT 

The instantaneous force acting on a circular cylinder 
immersed in a unidirectional periodic flovr characterized by 
U = V + sin 2irft was investigated. The total force was 
divided into three components consisting of a steady-state 
drag; an oscillating drag which is a function of the 
frequency parameter X = fd/V ; and an Inertial force which 
is also a function of the frequency parameter. The three 
resistance coefficients associated with the force components, 
Cj^(O), Cj^(X), and Cj^^(X) were calculated. 

It was found that the resistance coefficients are not 
materially affected by the free-stream oscillations whenever 
the frequency of the oscillation is lower than the natural 
frequency of the vortex shedding. However, when X approaches 
the Strouhal frequency of vortex shedding, the inertia 
coefficient is significantly reduced, and likewise when 
X approaches the value of twice the Strouhal number, Cp(X) 
is markedly reduced. In regions below this range, the 
Inertial force becomes negligible and the total resistance 
could be represented accurately by the coefficients C^CO) 
and Cj^(X) with fixed at its theoretical unseparated-flow 
value, 2.0. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

A. EXPERIMENTAL JUSTIFICATION 8 

B. SURVEY OF THE PREVIOUS INVESTIGATIONS 9 

C. SCOPE OP THE PRESENT STUDY I 5 

II. APPARATUS AND PROCEDURE I 7 

A. APPARATUS I 7 

1 . Water Tunnel and Oscillating Mechanism 17 

2. Test Specimens I 7 

3 . Instrumentation 21 

B. PROCEDURE 21 

1. Measurement of Force 21 

2. Measurement of Frequency 23 

3 . Measurement of Acceleration 23 

4. Measurement of Velocity 24 

5 . Tunnel Correction Factors 28 

III. THEORETICAL CONSIDERATIONS 32 

IV. EVALUATION AND PRESENTATION OF DATA 35 

A. RAW DATA 35 

B. EVALUATION OF COEFFICIENTS 40 

C. NORMALIZED PLOTS 43 

V. DISCUSSION OF RESULTS 50 

A. UNIDIRECTIONAL PERIODIC MOTION 50 

B. EXTENSION AND COMPARISON OP DATA WITH 

THOSE OBTAINED BY OTHERS 51 

VI. CONCLUSIONS 54 



3 



I 

I 

I 

I 



VII. RECOMMENDATIONS FOR FURTHER STUDIES 55 

LIST OF REFERENCES 56 

INITIAL DISTRIBUTION 'LIST 59 

FORM DD 1473 60 



4 



LIST OP FIGURES 



Figure 

1. Water Tunnel 18 

2. Butterfly Valve 

3 . Test Specimen 20 

4. Instrument Box and Cantilever Beam 22 

5 . Pressure Tap Diagram 25 

6 . Hydrodynamic Correction Factor 30 

7 . Theoretical Shape of the Squared Velocity 36 

8. Typical Chart Trace at Low Frequency of 

Oscillation 38 

9 . Typical Chart Trace at High Frequency of 

Oscillation 39 

10. Cp(0) for 10,000 < Re < 65,000 

11. Cp(0) for .0002 < X < .05 45 

12. Cp(X) for .0002 < X < . 05 46 

13 . Cp(0) for .01 < X < 2.5 r 4? 

14. Cp(X) for .01 < X < 2.5 48 

15 . C^(X) for .01 < X < 2.5 49 



5 



NOMENCLATURE 



Symbol 



A 


effective cylinder flow-blockage area (A = dA) 


a 


instantaneous acceleration 




cylinder cross-sectional area (A^^ = ird /4) 


Ats 


cross-sectional area of test section (A,pg = hJl) 


s 


steady state drag coefficient 


VO) 


coefficient for steady state portion of drag 


V^) 


coefficient for oscillating portion of drag 


Cj^(X) 


coefficient for inertia force 


d 


diameter of test cylinder 


e 


tunnel correction factor 


F 


instantaneous total force acting on the test 
cylinder 


f 


frequency of oscillation 


"d 


drag forces acting on the test cylinder 


L 


inertia forces acting on the test cylinder 


h 


height of- tunnel test section 


% 


constant (K^ = %pV^A) 




constant (K^ = W^Ap) 


K 3 


constant (K^ = Ap) 


h 


constant (K^ = PirfpA^ilV^) 




width of tunnel test section; test cylinder length 


n 


frequency of vortex shedding 


P 


pressure 


P 


height of dynamic pressure curve at t = T/4 
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q measured velocity 

R test cylinder radius 

r height of dynamic pressure curve at t = 3T/4 

Re Reynolds number (Re = Ud/v) 

Re Reynolds number (Re = Vd/v) 

St Strouhal number (St = nd/V) 

T period of oscillation 

t time 

U instantaneous velocity (U = V + sin2irft) 

V free stream velocity 

amplitude of velocity oscillation 

V average velocity 

X horizontal distance between pressure taps 1 and 3 

V horizontal distance between pressure tap 3 and 
center of test cylinder 

X normalized frequency of oscillation (X = fd/V) 

V fluid kinematic viscosity 

p fluid density 
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I . INTRODUCTION 



A. EXPERIMENTAL JUSTIFICATION 

The subject of forces acting on bluff bodies immersed 
in time-dependent flows has been and will continue to be of 
interest to fluid dynamicists, aerodynamicists , and practicing 
engineers for the special reason that in nature neither the 
body nor the fluid which surrounds the body is ever in a 
state of steady motion. Even the flow behind a bluff body 
moving steadily through a fluid is accompanied by large 
scale unsteadiness. Thus, any type of unsteadiness of the 
ambient flow and/or the motion of the body introduces 
additional changes in the characteristics of the flow and 
its analysis. 

The drag and inertial forces are interdependent as well 
as time dependent, and the resistance coefficients obtained 
in unseparated flows are not applicable to separated cases. 
Although indirect, the role of viscosity is paramount in 
that its consequences are separation, vortex formation and 
shedding, and resultant alternations in the virtual mass. 

It is thus clear that it is necessary to determine the 
relationships between various resistance components in terms 
of the unsteadiness of the ambient flow and/or body motion, 
geometry of the body, the degree of upstream turbulence, 
past history of the flow, etc. [Ref. 1]. 
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However, the understanding of the behavior of time- 
dependent flows about bluff bodies may be advanced only by 
considering the relatively more manageable cases and gradually 
Integrating the information so obtained. This led to the 
present experimental investigation of the unidirectional 
periodic flow about circular cylinders. 

B. SURVEY OF THE PREVIOUS INVESTIGATIONS^ 

A review of some of the previous investigations on time- 
dependent flov;s may best be presented by dividing the motions 
under consideration into various but admittedly arbitrary 
categories . 

1 . Theoretical Analysis of Unseparated Time-Dependent 

Laminar Flows 

The study of the theory of separation-free time- 
dependent laminar flows has enjoyed particular attention, 
partly due to its practical significance and partly due to 
its relative mathematical simplicity. Surveys by Stewartson 
[Ref. 2 ], Stuart [Ref. 3]> and Rott [Ref. 4] reflect the 
current level of understanding. 

The . unseparated class of unsteady flows, most of 
which results from the unidirectional or periodic acceleration 
of bodies in an Infinite or bounded fluid medium, gives rise 
to an induced mass which must be added to the real mass of 
the body. In general, the virtual inertia tensor depends 



^Thls section has been reproduced here, by permission, 
from an unpublished report by Prof. T. Sarpkaya. 
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on the shape of the body, the nature of the fluid medium, 
the orientation of the moving object with respect to the 
direction of motion, and the depth of relative submergence 
from a free surface and/or solid boundary. 

The effects of viscosity and separation on added 
mass or the dependence of added mass of a given body on the 
characteristics of the time-dependent separated flow cannot 
yet be theoretically evaluated nor can it be experimentally 
separated from the total resistance experienced by the body. 

2 . Small Amplitude Vibration of Bodies in a Liquid 
Otherv;lse at Rest and the Initial Instants of 
Acceleration of Bodies from a State of Uniform 
Velocity 

Vibratory motion has often been used [Refs. 5-8] as 
one of the experimental means to determine the virtual mass 
coefficients. The technique has been restricted to amplitudes 
of motion which are so small that either separation does not 
occur or the characteristics of separation do not change 
due to the imposed, high-frequency vibration. Thus, the 
results obtained are not applicable to occurrences in which 
the duration of unsteady flow in one direction is long enough 
for separation to occur. 

It is not necessary that the body be initially at 
rest to determine the added mass provided that one starts 
with a steady ambient flow about the body or with the body 
moving at constant velocity in a fluid othervjise at rest. 
Hamilton and Llndell [Ref. 9] have shown that the added-mass 
coefficient of a sphere determined by imposing a uniform 



10 



acceleration onto its steady state motion is equal to its 
theoretically determined value as well as to its experimen- 
tally determined value obtained by small-amplitude, high- 
frequency oscillations. 

Rayleigh [Ref. 10] has shown that if the acceleration 
continues, i.e. V = V + kt^, both the drag and the added-mass 
coefficients change with time. 

The foregoing arguments suggest that the force 
exerted by the fluid on a body immersed in a time-dependent 
flow might be expressed as 

F = ^ 2 CppA|v|V + dV/dt + history-dependent (1) 

drag and inertia forces 

in which A is the projected area of the body and m^ the mass 
of the fluid displaced by the body. The drag coefficient 
is to be taken equal to the quasi-steady state resistance 
coefficient and to be taken equal to that measured or 
calculated from the potential theory when the value of the 
history-dependent force is zero [Ref. 9l* 

3 . Bodies Subjected to Large Amplitude Harmonic 

Oscillations in a Fluid Otherv/ise at Rest 

A number of significant studies [Refs. 11-15] were 
made to determine the forces acting on bluff bodies under- 
going harmonic, large-amplitude oscillations in a liquid 
otherwise at rest. 
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Odar and Hamilton [Ref. 11] measured the force acting 

on a sphere subjected to harmonic oscillations and proposed 

a force equation comprised of three parts: a steady drag 

component where the drag coefficient is the well-established 

coefficient for steady translation of a sphere; an inertial 

component whose coefficient is determined experimentally as 

2 

a function of V /ad (a: acceleration; d: diameter of sphere); 
and a history-dependent component whose coefficient also is 
a function of V^/ad. 

H . Bodies Subjected to Wave Motion with Zero Mean Velocity 
Considerable work [Refs. l6-2l] has been done on this 
type of time-dependent flow in an attempt to predict the wave 
forces on piles and other submerged structures. In general, 
the total force acting on the body is assumed to be composed 
of a velocity-dependent and an acceleration-dependent force. 

The resulting equation, known as the Morison equation [Ref. 
l6], is given by 

F = JsCjjpAlVtV + dV/dt (2) 

where = 1 + C, C being the added mass coefficient. Experi- 
ments shov7 that Cp and depend on time and show considerable 

scatter about their assumed mean values [Refs. 17, 20] over 
a cycle. 
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5 . Bodies Subjected to Unidirectional Acceleration In 
a Fluid Otherwise at Rest or Unidirectional Unsteady 
Flow About Bodies Held at Rest 

Substantial effort has been made to determine the 
components of force acting on a body accelerating unldirec- 
tionally in a fluid otherwise at rest. Iverson and Balent 
[Ref. 22] and Keim [Ref. 23] towed spheres, disks, and cylin- 
ders in unsteady motion through still water. It was concluded 

that the total instantaneous force on these objects could be 

2 

described by a single coefficient of the form C = 2F/pAV 

and that this coefficient vras a function only of a so-called 

2 

acceleration modulus, ad/V . Keim, hov/ever, detected an 
apparent Reynolds number effect for cylinders. 

Laird et. al. [Ref. 24] measured the forces acting 
on accelerating cylinders and found strong evidence of devia- 
tion of the drag coefficient from the accepted values for 
uniform motion. They further found that the acceleration 
modulus did not correlate the resistance coefficient near 
boundary-layer transition. 

Sarpkaya and Garrison [Ref. 22] working with uni- 
directional flow with constant acceleration about cylinders 
and plates have found that both and depend on the 
history of motion, that the Morison equation (Eqn. 2) is 
valid only for flows with constant acceleration, and that 
the acceleration modulus can correlate the data only if the 
acceleration is kept constant. 
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6 . Forced or Self-Induced Transverse Oscillations of 
a Body in a Fluid In Steady Motion 

Experiments have been directed to the effects that 
attend the oscillations of a body in the plane of the lift 
force. This type of interaction becomes particularly impor- 
tant when the oscillation of the structure is the necessary 
condition for the generation of the exciting forces. The 
most pertinent information from these experiments is that 
the lift and drag forces act on the oscillating cylinder at 
the vortex shedding frequency and tv/lce the shedding fre- 
quency, respectively, provided that the driving frequency on 
the cylinder is appreciably different from the shedding fre- 
quency. V/hen the forcing frequency of the cylinder approaches 
the shedding frequency, the natural shedding frequency is 
lost and it "locks-in" to the forcing frequency. This 
synchronization persists over a range of frequencies which 
may be termed the "range of synchronization" [Ref. 26]. 

7 . Unidirectional Oscillatory Flow About Bodies at Rest 
and Steady Uniform Flovr About Bodies Subjected to 
Streamwlse Oscillations 

Relatively few studies have been conducted in this 
category. Chen and Ballengee [Ref. 27] examined the vortex 
shedding from circular cylinders in an oscillating free 
stream. Their results suggested that "in an oscillatory 
freestream of 3 HZ and Reynolds number up to 40,000, the 
vortex shedding from a circular cylinder responds instantan- 
eously to the freestream variations" and that "the instan- 
taneous Strouhal number stays sensibly constant at 0.2 ± 0.01." 
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Hatfield and Morkovin [Ref. 28] studied the effect 
of an oscillating free stream on the unsteady pressure on a 
circular cylinder. They have found that there is no signi- 
ficant coupling between the small-amplitude free-stream 
oscillations and the vortex shedding. Their results would 
suggest that the drag coefficient associated with the mean 
flow would essentially remain constant at its steady state 
value. On the other hand, Mercier [Ref. 29] who subjected 
cylinders to large streamwise oscillations found that the 
average drag coefficient significantly increases with fd/V 
and that the rate of Increase depends on the amplitude to 
diameter ratio. Combined, these results would suggest that 
the degree of coupling between the frequency of oscillation 
and the vortex-shedding frequency and hence the forces acting 
on the body strongly depend on the amplitude of oscillation. 

Davenport [Ref. 30] subjected bluff bodies to small 
amplitude oscillations in a water flume and evaluated the 
drag and inertia coefficients through the measurement of the 
rate of damping of the amplitude of oscillations of the 
bodies, l.e. v/ithout measuring the forces acting on the 
bodies. Using the frequency parameter fd/V to correlate 
the data, he found large variations in and C,j^. 

C. SCOPE OP THE PRESENT STUDY 

The present investigation is limited to an experimental 
determination of the instantaneous forces acting on a 
circular cylinder Immersed in a unidirectional periodic 
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flow characterized by 



U = V + sin 2TTft. (3) 

Its theoretical aspects are structured along the lines of 
the unpublished work by Sarpkaya [Ref. 313 conducted at 
Aerodynamische Versuchsanstalt-Gottingen v/ith respect to 
the supposition that the instantaneous force acting on an 
object immersed in a periodic flow strongly depends on the 
parameter fd/V and to varying degrees on W/v and 
where f is the frequency of oscillation and d the cylinder 
diameter. 

The range of these parameters during the experimental 
work v;as : 



0.0001 < fd/V <0.05 
11,000 < Vd/v < 63,000 
0.085 < V^/V <0.69 



As presented in Section IV, the parameter fd/V is extended 
to a value of 1.5 using the data of others. 
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II. APPARATUS AND PROCEDURE 

A. APPARATUS 

1. V/ater Tunnel and Oscillating Mechanism 

The experiments were conducted in a 500-gallon- 
capaclty, recirculating-type water tunnel as shovm in Figure 
1. V/ater was circulated by a low-rpm, high-capacity, centrif- 
ugal pump which was capable of providing test-section mean 
velocities of 2-25 feet per second. The test section 

I 

consisted of a rectangular conduit 4 inches wide, 8 inches 
high, and 16 inches long. 

Periodic variations in water velocity were provided 
by a rotating butterfly valve driven by a variable speed 
motor capable of a range of 0-i|0 rpm (Fig. 2). The magni- 
tude of the velocity fluctuations was dependent upon the 
loss in effective head through the valve. As this is partially 
a function of the valve-to-conduit area ratio [Refs. 35,36], 
an access plate was placed on the tunnel in order to change 
the effective valve area by adding concentric rings to the 
valve perimeter without removing or changing the valve body. 

Associated with the tunnel is a 15-gallon head tank 
located 12 feet above the test section centerline. 

2 . Test Specimens 

The test specimens consisted of polished aluminum 
cylinders placed across the width of the test section perpen- 
dicular to the direction of flow (Fig. 3) • A small bore 
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Figure 1: WATER TUNNEL 



ACCESS PLATE 

CONCENTRIC RINGS 

DIAMETER OP VALVE (including 
concentric rings) 

MOTOR AND REDUCTION GEAR 
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Figure 2: BUTTERFLY VALVE 
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Figure 3: TEST SPECIMEN 



hole drilled radially from the cylinder surface at the mid- 
section served as the total pressure tap. The test cylinder 
was mounted with one end held by a self-aligning bearing fit 
in one plexiglass window, with the other end passing through 
the other \-/indow to a cantilever beam force -measurement 
system. Two cylinders of diameters 0.75 and 1.0 inches were 
used. 

3 . Instrumentation 

The total force acting on the cylinder within the 
test section was measured directly by a strain-gage fitted 
cantilever beam (Fig. 4). The bridge circuit consisted of 
four active strain gages, with the output signal being fed 
to a dual-channel strip chart recorder. 

Three pressure taps were utilized for velocity and 
acceleration calculations. Two variable-sensitivity, 
differential-pressure transducers connected to two dual- 
channel strip chart recorders were used for this purpose. 

Each transducer/recorder system v;as calibrated with a water- 
manometer. 

B. PROCEDURE 

1 . Measurement of Force 

A continuous trace of the total, instantaneous force 
acting on the cylinder was obtained from the strain-gage 
fitted cantilever beam (Fig. which was wired to a strip 
chart recorder. This system was calibrated by static loading 
of the cylinder. Three different cantilever beams were used 
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FRONT VIEW 



B SELF-ALIGNING BEARINGS 
CB CANTILEVER BEAM 

IB INSTRUMENT BOX 

SG STRAIN-GAGES 
TS TEST SECTION 



Figure 4: INSTRUMENT BOX AND CANTILEVER BEAM 



during the course of the investigation, each beam being 
several times more sensitive than the previous one, in order 
to increase the range of experimentation. 

2 . Measurement of Frequencies 

Each strip chart recorder has the capability of 
variable, controlled paper speeds. Thus, with chart speed 
and distance betvjeen corresponding points on the trace 
kno^'m, the period of the cyclic signal could be determined. 
The frequency of oscillation was then simply the reciprocal 
of the period. 

3 . Measurement of Acceleration 

Acceleration, with the exception of a scale factor, 
was measured directly from the differential pressure taken 
from two static pressure taps placed a knovm distance apart. 
Neglecting friction losses, the equation of motion yields 
(Fig. 5), 



(P-j_ - P.^)Arpg = ma 

where m is the mass of fluid, 
one has 



( 4 ) 

Writing m = p X A^^ in (4), 



(P^ - P 3 )A ^2 = p X A^g a 



(5) 



or 



a = AP^_^/pX. 



( 6 ) 
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4 . Measurement of Velocity 



2 

The dynamic pressure, in the form pU /2, was 
determined using the differential pressure from the total 
and static pressure taps 2 and 3> respectively, (Fig. 5), 
by properly locating the wall pressure tap for a given size 
cylinder. 

Determination of this distance begins with the 
unsteady Bernoulli equation 

DV ^ ^ = _ 1 iL (7\ 

Dt8t''8s~p8s ^ 

which reduces to 



(P + pvVs) = P ||. 



( 8 ) 



Upon integration, (8) becomes 



(P + pV^/2) 



= -P / ^ ds 
1 1 



(9) 



or 



? ? ^ Pi V 

?2 + PV 2 /2 - - PV^ /2 = -p / ds. 



(10) 



It is evident from Figure 5 that V 2 = 0 at the stagnation 
point. Therefore, (10) becomes 



^^1-2 " ^2 " ^1 " "p If p^^/2. 



2H 



( 11 ) 
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Figure 5: PRESSURE TAP DIAGRAM 



Now, from Equation ( 6 ) 



' ’’3 - ''1 ' H 



( 12 ) 



and combining these pressure differences yields, 



P 2 - P 3 = pvVp - p^/ ^ ds + px II . (13) 

2 

In order to satisfy the requirement that P^ - P 3 = pV /2, 
it is necessary, considering ( 13 ), that 



X 



9t 




ds . 



(14) 



To evaluate this integral, the variation of V with s along 
the stagnation streamline must be known. 

The complex potential for a cylinder in a uniform 
stream is given by 



= V(Z + 




(15) 



and the complex velocity is 



W 



dZ 




But, by definition 



( 16 ) 



VJ = u - iv . 



(17) 
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Therefore, along the stagnation streamline connecting the 
points 1 and 2, one has 



V = 0 
u = V(1 




(18) 



since Z = -s. Substituting this result into the integral 
of Eqn. (l4) yields 



/ 

1 




ds 



at 



f il - ^) ds. 
1 s 



(19) 



Carrying out the integration and equating the results with 
the requirement of (l4) yields 



X 



at 



av, ^ R^x2 

— ^1 



( 20 ) 



which simplifies to. 



X 



(X + Y) ,+ 



R 

(X + Y) 




( 21 ) 



or 

^ TX^) ' 2R (22) 

where (X + Y) is the distance to the point 2 and R is the 
radius of the cylinder and distance to the point 1. 
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Since the flow upstream from the cylinder extends 



to relatively large distances, i.e., as X ^ Equation (22) 
reduces to 



Y = 2R. 



(23) 



Therefore, pressure tap 3 (Pig* 5) must be located one 
diameter upstream from the center of the cylinder in order 
to obtain P2 - PlT“/2. 

The length of the test section upstream from the 
cylinder was sufficiently large to Ignore the small influence 
of the contraction of the flow near the test section on the 
ambient flow. It is for this reason that X is taken to be 
effectively Infinitely large. 

5 . Hydrodynamic Tunnel Correction Factors 

Steady state tunnel correction factors as discussed 
in Pope [Ref. 3^] account for the fact that the bluff body 
is not being tested in an infinite medium. There is a 
blockage-effect correction to the measured velocity, which 
in turn affects the calculated drag coefficient. In a usable 
form, this correction factor is given by 



e 



I2 ^h^ ^ 



(24) 



where d is the cylinder diameter and h is the height of the 
tunnel test section. This correction factor is then used to 
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modify the measured velocity, q, by 



U = q(l + e) 



( 25 ) 



which then modifies the uncorrected drag coefficient to 



= 2P/pU^£d . 

^corr 



( 26 ) 



The effect of this correction is to increase the measured 
velocity and thus to decrease the drag coefficient. 

However, with the critical placement of the static 
pressure tap 3 and its proximity to the test cylinder, 
another correction must be applied to the measured velocity. 
Prom Pigure 6, for steady state flows, the ratio of local 
fluid velocity to the free stream velocity is given by 



U = 



2 If 

1 — 2 — ^ cos 20 + 

r r 



(27) 



Por the geometry in this situation, the effect of this 
correction is to decrease the measured velocity and thus 
to increase the drag coefficient. 

As the effects of (25) and (27) counteract each other 
in the steady state flov; situation, and since no information 
is available for time-dependent flows, neither correction 
was applied. This is particularly justified in the present 
study since the time-dependent drag coefficients were compared 
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Figure 6: HYDRODYNAMIC CORRECTION FACTOR 
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with calculated steady drag coefficients, neither being 
corrected for any hydrodynamic constraints. 
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III. THEORETICAL CONSIDERATIONS 



Assuming that the unidirectional periodic flow in which 
the test cylinder is immersed is a sinusoidal flow superposed 
on a steady flow, the instantaneous free stream velocity can 
be represented as 



U = V + sin 2 it ft. 



( 28 ) 



Therefore, the instantaneous acceleration becomes 



a = ^ = 2TTfV. cos 2 tt ft 
at A 



(29) 



and the velocity squared is 



+ 2W^ sin 2tt ft + sin^ 27rft. (30) 



Now, in dealing with the total force acting on the 
cylinder in time-dependent flows, it has been most advan- 
tageous to consider it as composed of a drag force and an 
inertial force, such as 



F = Fp + Fj. (31) 

Then using the concept of the normalized drag and inertia 
resistance coefficients, C^ and C^^, respectively, the total 
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force can be represented as 



P = 0p%Apu2 + § (32) 

where A Is the effective cylinder flow-blockage area, 
is the cylinder cross-sectional area, I is the cylinder 
length, and p the fluid density. Upon substitution of (29) 
and ( 30 ) into (32), one obtains, 

P = Cp[%ApV^ + ApW^sin2iTft + %ApV^^sin^2Trft ] 

(33) 

+ Cj^^[p£A^2TTfV^cos2TTft]. 

However, the results of the works previously cited indicate 
that a two parameter total-force equation does not seem to 
be capable of providing sufficient generality to cover the 
conditions encountered in arbitrary time-dependent flows. 
Also, the past history of the flow must be considered along 
with the velocity, rate of change of velocity, and Reynolds 
number. Therefore, defining Cj^(O) as the steady state 
coefficient of drag v/hich occurs at time t = 0 when U = V, 
and assuming that a normalized frequency parameter X = fd/V 
provides the necessary correlation between the resistance 
coefficients and the aforementioned flow parameters, the 
total force can be written as 
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F = Cp(0)J5ApV^ 

+ Cj^(X)[W^Apsin2Tift + %ApV^^sin^2Trft ] (3^) 

+ Cj^( X ) [A^p-C2T7fV^cos2Tift ] . 

The coefficients Cj^(O), Cj^(X), and Cj^(X) are determined 
experimentally. It is obvious that the validity of the 
foregoing analysis depends upon the existance of a correla- 
tion between the said coefficients and the normalized 
oscillation parameter, X, for various values of Reynolds 
number, X, and 
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IV. EVALUATION AND PRESENTATION OP DATA 



A. RAW DATA 

As cited earlier, the differential pressure directly 
provides the velocity in squared form. Figure 7 shows a 
theoretical curve of = (V + V^ sinPnft)^ with V = 1.5 
and V^ = 0 . 5 . The distances labeled p and r occur at 
t = T/^ and t = 3T/^, respectively. With these distances 
known, V, V^, and U at time t = T/8 can be calculated 
as follows. Writing at the specific times. 



one has 



and 



t = t/4 


= (V + 


Va)= = P 


(35) 


t = 3T/4 


= (V - 


u 

11 

CM 

<c 

> 


(36) 


t = T/8 


= (V + 


V^ \l2/2)‘^ 


(37) 



V + V^ = sj^ 


(38) 


Cl 

1 

c 

> 

11 


(39) 



which yields 



V = ( \/p + sJt)/2. 



(40) 
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Figure 7: THEORETICAL SHAPE OF THE SQUARED-VELOCITY 



Now, upon subtracting (36) from (35), one obtains 



p - r = (V + V^)^ - (V - V^)^ = i}W^ (ill) 

which reduces to 

= (p - r)/W. (i»2) 

Also upon expansion of equation (37), at time t = T/8, one 
has 



+ \y^VV^ + (43) 

Figures 8 and 9 display copies of strip chart recordings 
of typical data sets for low and high frequencies of oscilla- 
tion, respectively. The upper traces represent 
2 

AP 2_2 “ pU /2 vjhlle the lower traces Indicate the Instantan- 
eous total force.. Due to the fact that the torque curve of 
a butterfly valve Is not symmetrical about 45° , nor the 
discharge coefficient (and thus the head loss) symmetrical 
about 90° In unidirectional flow [Ref. 33], the velocity 
curve Is not a pure sinusoid. Therefore, the following 
procedure was used In reducing the raw data to usable form: 

1. The period, T, was measured between any two corres- 
ponding points on the curve. Thus T was the true period of 
the repetitive signal. 
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Figure 8: TYPICAL CHART TRACE AT LOW FREQUENCY OF OSCILLATION 
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Figure 9: TYPICAL CHART TRACE AT HIGH FREQUENCY OF OSCILLATION 



2. The distances p and r were measured to the highest 
and lowest points on the pressure curve, respectively. 

The occurrance of these points were labeled T/^ and 3T/4 
Irrespective of their distance on the time axis since the 
minimum and maximum velocities will occur at these times. 

3. The height at the squared-velocity curve at times 
t = 0 and t = T/8 were calculated from, p and r as If the 
velocity were a true sinusoid. The corresponding points on 
the curve were then found at these heights, and the base 
line marked 0 and T/8. 

. Vertical lines were then drawn from the pressure 
curve to the force curve at the places marked 0, T/8, and 
T/i|, and the height of the force curve tabulated at these 
points. Since the two traces are made simultaneously and 
represent Instantaneous values, these quantities then give 
the forces which would have occurred at times 0, T/8, and 
T/il had the velocity been a true sinusoid. 

5 . The values of velocity and force were then calcu- 
lated from the calibration equations. Reynolds number 
based on the mean velocity v;as also calculated. For this 
purpose, an average water temperature of 70° P was assumed 
for evaluating the density and viscosity of the tap water 
In the tunnel. 

B. EVALUATION OP COEFFICIENTS 

As detailed In Part III, the total force acting on the 
cylinder Is assumed to take the form 



Ho 



F = Cp(0)[%Apv2] 



+ Cp(X)CW^Apsin2TTft + ApV^^sin^2Tift ] {HH) 



+ Cj^(X)[2TTA^p£fV^cos2Trft] 



After the reduction of the raw data, at times t = 0, t = T/8, 
and t = T/4 all the variables with the exception of Cj^(O), 
Cj^(X), and Cj^(X) are known. However, a system of three 
equations for the three unknowns can be developed. Re- 
writing (4i|) as 



and evaluating (^5) at times 0, T/4, and T/8 results in 



F = Cp(0)K^ 



+ Cp(X)[K2Sin27Tft + K2Sin^2Trft] 



(^5) 



+ Cj^^(X) [K|^cos2irft ] 



t = 0 F = Cp(0)K^ + Cj^(X)K^ 



(il6) 



t = t/ 8 F = Cjj(0)K^ + Cjj(X)[K2 \/2/2 + (^7) 

+ Cj^(X)Ki^ J 2/2 



t = T/ll F = Cj^(0)K^ + Cp(X)[K2 + K^] 



(il8) 
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These three equations can be solved simultaneously for Cp(0), 
Cp(X), and Cj^(X). 

However, in the range of experimentation of this 
study, the ratio of inertia forces to that of the drag 
forces was << 1. This fact led to problems of instability 
in the matrix used to solve for the inertia coefficient. 

For this reason, Cj^(X) was fixed at a chosen value, and the 
values of Cp(0) and Cp(X) were solved using the instantaneous 
data at times t = 0 and t = T/^. The two resulting equations 
were 



t = 0 



Cp(0) 



p - 







(iig) 



and 



F - C^(0)K^ 

The coefficients 0^(0) and Cp(X) were calculated from 
(^9) and (50) by using two different values for Cj^(X). 

First, a theoretical value of Cj^(^) = 2 for unseparated 
potential flow and second, a semi-theoretical and semi- 
experimental value of Cjyj(X) = 1.5 [Refs. 25,35] were used. 
The resulting values of 0^(0) and Cj^(X) were essentially 
unchanged. The difference between Cj^(O) with Cj^(X) = 2 
and Cj^(O) v;ith Cj^(X) = 1.5 was negligible to the second 
decimal. This also held true for Cp(X). 
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C. NORMALIZED PLOTS 

The non-dimensional drag coefficient Cp(0) is pre- 
sented in Figure 10 as a function of the Reynolds number 
based on average velocities. 

The coefficients Cj^(O), Cj^(X), and Cj^(X) are then 
presented as functions of the normalized frequency parameter, 
X (Pigs. 11-15). 
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Figure 10: C^(0) VERSUS REYNOLDS NUMBER 
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Figure 12: C„(X) FOR .0002 
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Figure 14: 
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V. DISCUSSION OF RESULTS 



A. UNIDIRECTIONAL PERIODIC MOTION 

Although the flow remained unidirectional and periodic, 
its oscillatory motion was not quite sinusoidal, as discussed 
in Section IV. This aberration in the flov: appeared to 
manifest itself differently for high and low frequencies of 
oscillation. Irrespective of the value of X. 

At low frequencies of oscillation (Fig. 8) the curve of 

2 

pU /2 appears symmetrical about the vertical with the abso- 
lute values of acceleration and deceleration about equal. 
However, there is a long period in the cycle at which the 
valve appears to have no effect on the flow. At the highest 
portion of the curve the valve is in the horizontal position 
with respect to the conduit axis (minimum head loss), while 
at the lowest portion of the curve the valve is vertical 
(maximum head loss). It is thus apparent that the available 
head is not significantly changed when the valve is about 
± 30 degrees from the horizontal, causing a flat portion 
in the curve. 

At high frequencies of oscillation (Fig. 9), this ’’null 
portion" of valve travel creates a pressure curve which is 
unsymmetrical about the vertical characterized by a small 
acceleration and rapid deceleration. 

Several attempts were made to alleviate this problem. 
First, a baffle plate was installed upstream of the valve 
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to decrease V. At this time vertical ridges extending 
downstream were attached around the holes in the baffle in 
order to decrease the magnitude of the deceleration. Second, 
small vertical plates were attached perpendicular to the 
valve faces so that the valve might affect the head loss in 
the near horizontal position. Finally, two plates were 
attached perpendicular to the valve faces with the same 
diameter as the valve itself to form a four-bladed butterfly 
valve. Unfortunately, time did not permit manufacturing the 
plates to the same width or contour as the existing valve. 

It is for this reason that only those portions of a 
cycle which most nearly exhibited a sinusoidal velocity 
variation ivere used in evaluating the various coefficients. 

A butterfly valve system consisting of several plates is 
anticipated to alleviate this problem and yield a fairly 
sinusoidal oscillation in velocity over the entire cycle. 

B. EXTENSION AND COMPARISON OP DATA WITH THOSE OBTAINED 

BY OTHERS 

The data obtained herein are presented together with 
those obtained by Davenport [Ref. 36 ] at the University of 
Bristol and by Sarpkaya [Ref. 31] at the Aerodynamlsche 
Versuchsanstalt , Gottingen, using a circular cylinder of 
1 inch diameter in an oscillating wind tunnel within the 
range of Reynolds numbers from 5000 to 90,000. 

First, it is apparent that the data of Davenport and 
Sarpkaya constitute extensions of the present data to higher 
values of the frequency parameter, X. 
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Second, It is clear that all of the data correlate fairly 
well with X over a large range of X and Reynolds numbers 
regardless of whether it is obtained with air or water in 
a small or large tunnel. 

Third, it is evident that the effects of the free-stream 
oscillations on Cj^(X) become significant only when the 
frequency of oscillation approaches the natural vortex- 
shedding frequency, l.e. when X approaches the Strouhal 
number of about 0.2. Although a theoretical explanation of 
the decrease of Cjyj(X) as X -»■ St is beyond the scope of the 
present study, it is not too much of a conjecture to state 
that the formation, growth, and motion of the vortices are 
significantly altered by the oscillations when the shedding 
frequency is in the neighborhood of the oscillatory frequency 
of the free stream. \Vhen the period of oscillation is several 
orders of magnitude larger than the period of vortex shedding, 
the flow behaves essentially as a juxtaposition of instan- 
taneous quasi-steady states and the total resistance can 
be represented accurately by the coefficients Cj^(O) and 
Cp(X). 

Finally, it is evident that Cj^(X) decreases sharply in 
the range of fd/V from 0.25 and 0.5. The reasons for this 
are as follows. It is a well-known fact that the drag 
acting on a cylinder undergoes oscillations with a frequency 
twice the frequency of vortex shedding. In other words, as 
the lift force changes its direction, up or dov;n, the drag 
force oscillates with an amplitude of AC^ about a mean value 
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of C_. Denoting the oscillation frequency of drag with 
fp, we then have fj^d/V ~ 0.44. The same parameter for the 
lift oscillations is ~ 0.22 vjithin the range of 

Reynolds numbers studied. Thus the data presented in Fig. l4 
show that Cjj(X) is, and should be, affected in the range 
of fd/V from 0.25 to 0.5j i.e. as the drag-oscillation 
frequency nearly synchronizes v;ith the free-stream oscillation 
frequency. Only a complete analytical investigation, if 
at all possible, can explain as to why Cp(X) should decrease. 
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VI . CONCLUSIONS 



The conclusions which can be drawn from this study are : 

1. The frequency parameter X correlates fairly well the 
data obtained by various investigators regarding the resis- 
tance coefficients. 

2. In the range of 0 < X < 0.05> the resistance 
coefficients are not materially affected by the free-stream 
oscillations. This range has not been previously investi- 
gated. All or most of the previous studies have dealt with 
oscillations v;ith a zero mean velocity, i.e. with X ->■ «>. 

3. The drag and inertia coefficients vjhich are functions 
of the frequency parameter X are significantly reduced in 
the range where the free stream oscillations nearly coincide 
with the flow induced oscillations of lift and drag. 

4. The reasons leading to the decrease in Cj^(X) and 
Cj(^(X) in the vicinity of X = 0.4 and X = 0.2 respectively, 
can be clarified only by a combination of additional theoret- 
ical and experimental investigations where the full Navier- 
Stokes equations are solved for an oscillatory flow about a 
circular cylinder and where the characteristics of flow in 
the near-wake of the cylinder are measured. 
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VII. RECOMMENDATIONS FOR FURTHER STUDIES 



The following are recommendations for further studies 

1. Improve the butterfly valve system so that nearly 
sinusoidal oscillations may be obtained in the range of 

0 < X < 2.0. 

2. Undertake a study of the solution of the Navier- 
Stokes equations via finite-difference or finite-element 
methods for the range of X values cited above. 

3. Carry out experiments with flat plates and other 
shapes of bodies to ascertain the universality of the 
parameter X. 
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